function testlomba()

% plot results of lomb amplitude (denormalized) versus input AMPLITUDE
% across a few frequencies.

h = findobj('Tag', 'testspecx');
if(isempty(h))
   h = figure('Tag', 'testspecx', 'name', 'test spectral analysis', 'NumberTitle', 'off');
end;
figure(h);
clf;
freq = 100;
amp = 1;
flomb = [1:300];

freq = [20,30,50,60,70,100];
a = [0.1,0.2,0.5,1,1.5,2,3];
dur = [25, 30, 50, 75, 100, 200];
for i = 1:length(a)
   [fl(i,:), pvl(i,:), Y] = testspec2(a(i),0,freq(3),dur(4));
   [mx(i), pf(i)] = max(pvl(i,:));
end;
figure(h);
plot(a, mx, 'ro');
b = mx./a
c = flomb(pf)./freq


